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Abstract 

We present results from a 150 ksec Suzaku observation of the Seyfert 1.5 NGC 
3516 in October 2005. The source was in a relatively highly absorbed state. Our 
best-fit model is consistent with the presence of a low-ionization absorber which has 
a column density near 5 x 10^^ cm~^ and covers most of the X-ray continuum source 
(covering fraction 96-100%). A high-ionization absorbing component, which yields a 
narrow absorption feature consistent with Fe K XXVI, is confirmed. A relativistically 
broadened Fe Ka line is required in all fits, even after the complex absorption is taken 
into account; an additional partial-covering component is an inadequate substitute 
for the continuum curvature associated with the broad Fe line. A narrow Fe Ka 



emission line has a velocity width consistent with the Broad Line Region. The low- 
ionization absorber may be responsible for producing the narrow Fe Ka line, though 
a contribution from additional material out of the line of sight is possible. We include 
in our model soft band emission lines from He- and H-like ions of N, O, Ne and 
Mg, consistent with photo-ionization, though a small contribution from coUisionally- 
ionized emission is possible. 

Key words: galaxies: active — X-rays: galaxies — galaxies: individual (NGC 
3516) 

1. Introduction 

The 6.4 keV Fe Ka emission line has long been known to be an important diagnostic of 
the material accreting onto supermassive black holes. The Compton reflection hump, frequently 
seen in Seyfert spectra above ~7 keV and peaking near 20-30 keV (Pounds et al. 1990), indicates 
that Seyferts' Fe K lines may have an origin in optically thick material, such as the accretion 
disk. Observations with ASCA indicated that many Fe Ka lines were broad (FWHM velocities 
up to ~0.3c) and asymmetrically skewed towards lower energies, implying an origin in the inner 
accretion disk; the line profile is sculpted by gravitational redshifting and relativistic Doppler 
effects (e.g., Tanaka et al. 1995, Fabian et al. 2002). However, XMM-Newton and Chandra 
observations have been revealing a more complex picture. A narrow Fe K component (FWHM 
velocities ~5000 km s~^ or less) appears to be quite common; these lines' widths suggest 
emission from distant material, such as the outer accretion disk, the optical/UV Broad Line 
Region (BLR) or the molecular torus. Spectral observations in which the broad and narrow 
components are deconvolved are thus a prerequisite for using the Fe K line as a tracer of the 
geometry of the emitting gas. 

At the same time, there is strong evidence from X-ray and UV grating observations for 
the presence of ionized material in the inner regions of a large fraction of AGN (e.g., Blustin 
et al. 2005; McKernan et al. 2007). High-resolution spectroscopy shows the gas is usually 
outflowing from the nucleus; typical velocities are ~ a few hundred km s~^. Absorption due 
to a broad range of ionic species is commonly seen; and for many sources, there is evidence for 
several different photo-ionized absorbing components, as opposed to a single absorber, along 
the line of sight. In the Fe K bandpass, some Seyferts show evidence for absorption by H- or 
He-like Fe, indicating a zone of highly-ionized absorbing material (e.g., NGC 3783, Reeves et 
al. 2004). 

Cold absorbing gas, with line of sight columns in excess of the Galactic value, are rou- 
tinely observed in Seyfert 2 AGN in accordance with unification schemes (Urry, Padovani 1995), 
and have also been reported in some Seyfert 1 AGN. Importantly, variations in column density 
on timescales from hours to years have been observed in both Seyfert 1 AGN (e.g., in I Zw 1, 
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Gallo et al. 2007; see also Lamer et al. 2003, Puccetti et al. 2004) and Seyfert 2 AGN (Risaliti 
et al. 2002; Risaliti et al. 2005), suggesting that the absorbing circumnuclear material is not 
homogeneous in either Seyfert type, has a high tranverse velocity and occurs over a range of 
length scales. 

NGC 3516 is a well-studied, nearby (z = 0.008836; Keel 1996) Seyfert 1.5 AGN that 
can display strong 2-10 keV flux variability on timescales of hours to years (e.g., Markowitz, 
Edelson 2004). Previous X-ray spectroscopic observations of NGC 3516 e.g., Nandra et al. 
(1997) using ASCA, have indicated the presence of a broad Fe K line, and this source is known 
to also contain complex and ionized absorption. Numerous UV absorption lines, including N 
V, C IV and Si IV, were observed with the International Ultraviolet Explorer (Ulrich, Boisson 
1983); absorption line strengths vary on timescales as short as weeks as the absorber responds 
to variations in the ionizing flux (e.g., Voit et al. 1987). Hubble Space Telescope observations 
have revealed that this component of absorbing gas (henceforth called the "UV absorber" ) may 
consist of several distinct kinematic components (Crenshaw et al. 1998). 

X-ray spectra of NGC 3516 can exhibit evidence for large columns (^10^^ cm~^) of 
absorbing gas (e.g., Kolman et al. 1993), and the X-ray absorbers can display variability on 
timescales of years (e.g., Mathur et al. 1997). Using Chandra gratings data from observations 
in April 2001 and November 2001, Turner et al. (2005) observed K-shell absorption lines due 
to H- like Mg, Si and S, and He-like Si, evidence for a highly-ionized absorber, likely with 
column density ^10^^ cm~^, outflowing at ~1100 km s^^. Simultaneous with these Chandra 
observations in 2001 were two XMM-Newton observations. Turner et al. (2005) modeled the 
continuum curvature of the two XMM-Newton EPIC spectra by including a partial covering, 
mildly- ionized absorber; the column density was ~2.5 x 10^^ cm~^, with a covering fraction 
of ~50%. However, the formal requirement for the broad Fe line was reduced, leading to 
uncertainty as to whether the broad Fe line really existed in NGC 3516. Spectral fitting using 
an instrument with a wide bandpass is thus necessary to remove such model degeneracies. 

In this paper, we report on an observation of NGC 3516 made with the Suzaku obser- 
vatory in October 2005. The combination of the X-ray Imaging Spectrometer (XIS) CCD and 
the Hard X-ray Detector (HXD) instruments have yielded a broadband spectrum covering 0.3 
to 76 keV, allowing us to deconvolve the various broadband emitting and absorbing compo- 
nents. Furthermore, the exceptional response of the XIS CCD and high signal-to-noise ratio 
of this observation have allowed us to study narrow emission lines in great detail. §2 gives a 
brief overview of the Suzaku observatory, and describes the observation and data reduction. 
§3 describes fits to the time-averaged spectrum. Variability analysis is briefly discussed in §4. 
Flux-resolved spectral fits are discussed in §5. In §6, we describe a search for narrow red- or 
blue-shifted lines in the Fe K bandpass. The results are discussed in §7, and a brief summary 
is given in §8. 
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2. Observations and Data Reduction 



The nucleus of NGC 3516 was observed by Suzaku from 2005 October 12 at 13:45 UT 
until October 15 at 09:07 UT. Suzaku was launched 2005 July 10 into a low-Earth orbit. It 
has four X-ray telescopes (XRTs; Serlemitsos et al. 2007), each with a spatial resolution of 2' 
(HPD). The XRTs focus X-rays onto four X-ray Imaging Spectrometer (XIS; Koyama et al. 
2007) CCDs, which are sensitive to 0.2-12 keV X-rays on a 18' by 18' field of view, contain 
1024 by 1024 pixel rows each, and feature an energy resolution of ~150 eV at 6 keV. Three 
CCDs (XISO, 2 and 3) are front-illuminated (FI), the fourth (XISl) is back-illuminated (BI) 
and features an enhanced soft X-ray response. The XRT/XIS combination yields effective areas 
per detector of roughly 330 cm^ (FI) or 370 cm^ (BI) at 1.5 keV, and 160 cm^ (FI) or 110 cm^ 
(BI) at 8 keV. Each XIS is equipped with two ^^Fe calibration sources which produce fluorescent 
Mn Ka and K/3 lines and are located at the CCD corners. Suzaku also features a non-imaging, 
collimated Hard X-ray Detector (HXD; Takahashi et al. 2007); its two detectors, PIN and GSO, 
combine to yield sensitivity from ~10 to ~700 keV. Further details of the Suzaku observatory 
are given in Mitsuda et al. (2007). 

2.1. XIS Reduction 

The XIS data used in this paper were version 1.2 of the screened data (Fujimoto et al. 
2007) provided by the Suzaku team. The screening is based on the following criteria: grade 
0, 2, 3, 4, and 6 events were used, the cleansis script was used to remove hot or flickering 
pixels, data collected within 256 s of passage through the South Atlantic Anomaly (SAA) were 
discarded, and data were selected to be 5° in elevation above the Earth rim (20° above the 
day-Earth rim). The XIS-FI CCDs were in 3x3 and 5x5 editmodes, for a net exposure time 
after screening of 135.0 (XISO), 134.8 (XIS2) and 135.2 (XIS3) ks. XISl was also in 3x3 and 
5x5 editmodes, for a net exposure of 135.4 ks. The XIS was in normal clocking mode. 

The source was observed at the nominal center position of the XIS. For each XIS, we 
extracted a 3' radius centered on the source. The background was extracted using four circles 
of radius 1.5', each located ~6' from the source, but chosen to avoid the z = 2.1 QSO RX 
J110741.4-F 723235, located 4.5' SE of NGC 3516. Spectra were binned to a minimum of 25 
counts bin~^ to allow use of the statistic. 

Response matrices and ancillary response files (ARFs) were generated for each XIS in- 
dependently using XISSIMRMFGEN and xissimarfgen version 2006-10-26 (Ishisaki et al. 2007). 
The ARE generator takes into account the level of hydrocarbon contamination on the optical 
blocking filter. We estimate a carbon column density of 0.8, 1.2, 1.7, and 2.8 x 10^® cm~^ for 
XISO, 1, 2 and 3, respectively. Finally, we co-added the three XIS-FI spectra using mathpha, 
and co-added the response files and ARFs using addrmf and addarf, respectively. 

To examine the accuracy of the XIS RMFs and determine residual line width due e.g., 
to imperfect CTI correction, we generated spectra for the ^^Fe calibration source lines on each 

4 



XIS using the above response matrices and ARFs. We fit the cahbration source spectra with 
three Gaussians. Two Gaussians were for the Mn Ka doublet (expected energies 5.899 keV 
and 5.888 keV), with energy centroids fixed to be 11 eV apart, and the higher energy hne flux 
set to twice that of the lower energy one. The third Gaussian was used to model the K/? line, 
expected at 6.490 keV. We found the average of all the calibration line widths cr to be < 4 
eV. The Mn Ka line energy centroids for the co-added FI spectrum were consistent with the 
expected energies to within 1 eV. For XISl, the Mn Ka line energy centroids were 3±2 eV 
lower than expected. Such discrepancies are well within the accuracy (~0.2% at Mn-Ka) of 
the energy calibration of XIS. Fitting the calibration source lines without the response file, we 
determined the FWHM energy resolution during the observation to be 145 eV (average of the 
4 XlSes). 

2.2. HXD Reduction 

We used data from the HXD-PIN only; NGC 3516 was extremely faint in the HXD-GSO 
band, and some aspects of the GSO background are still being studied, so we defer analysis of 
the GSO data to a later time. The PIN source spectra were extracted from cleaned version 1.2 
(pre-1.2-rl) HXD event files provided by the HXD instrument team. PIN background count 
rates are variable and strongly depend on the time since SAA passage (Kokubun et al. 2007), 
so we selected data according to the following criteria: at least 500 s since SAA passage, cutoff 
rigidity (COR) > 8 GV, and day- and night-Earth elevation angles each >5°. Instrumental 
(non- X-ray) background spectra for the PIN were provided by the HXD Team ("Background 
A" model) generated from a time- dependent model. The model utilized the count rate of upper 
discriminators as the measure of cosmic ray flux that passed through the silicon PIN diode and 
yielded background spectra based on a database of non X-ray background observations with the 
PIN (Fukazawa et al. 2007). The current accuracy of the PIN non-X-ray background (NXB) 
model for a 1 day observation is about 5% (peak-to-peak residuals). Both the source and 
background spectra were generated with identical good time intervals, and the source exposure 
was corrected for instrument dead time (a ~5% effect). This yielded a good time exposure of 
105.8 ks. 

Data < 12 keV were discarded due to noise contamination near the lower threshold 
of the PIN diode. Data above 76 keV were also discarded: the gain above an internal Bi Ka 
calibration line at 76 keV is not well-defined, though the photon statistics above this energy were 
poor anyway for this observation. Further details of the HXD in-orbit performance are given 
in Kokubun et al. (2007). To model the contribution to the total background from the Cosmic 
X-ray Background (CXB), a spectrum of the form 9.0 x 10'^{E/3keV)-°-^^ exp(-^/40keV) 
erg cm~^ s~^ sr~^ keV~^ (Gruber et al. 1999) was used. We note that some recent works (e.g., 
Frontera et al. 2007) suggest a 10% normalization increase compared to Gruber et al. (1999). In 
addition, spatial fluctuations of order ~5-10% over ~ a square degree are known (e.g., Barcons 
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et al. 1998). However, the effect on the net spectrum was neghgible; for instance, the change 
in Compton reflection component strength was 1%, typically. To simulate a CXB spectrum 
using XSPEC, we assumed a model of the form powerlaw*highecut^, with photon index T = 
1.29, cutoff energy of 0.1 keV, and an e- folding energy of 40 keV. The power-law normalization 
of 8.8 X 10~^ ph cm~^ s~^ keV~^ (at 1 keV) was used, appropriate for a source observed in 
XlS-nominal mode. The 12-76 keV CXB flux was 1.4 x 10"^^ erg cm~^ s~^ (using the Gruber 
et al. 1999 normalization). The total (X-ray plus particle) background 12-76 keV flux was 4.4 
X 10^^° erg cm"^ s~^. 

The spectrum was binned to a minimum of 400 count bin~^. We used the response file 
ae_hxd_pinxinom_20060814.rsp. The mean 12-76 keV net source flux and count rate were 1.1 
X 10~^° erg cm~^ s~^ and 0.16 ct s~^, respectively. Figure 1 shows the net source, background, 
and total (source + background) spectra. The source spectrum is always at least 15% of the 
total up to ~40 keV. 

3. Model Fits to the Time- Averaged Spectrum 

We used 0.4-11.5 keV data in the XIS-FI spectrum and 0.3-10.0 keV data in the XIS-BI 
spectrum. We ignored 1.72-1.87 keV in the co-added FI spectrum and 1.80-1.87 keV in the 
BI spectrum due to uncertainties in calibration associated with the instrumental Si K edge. 
In all fits, we included a constant to account for relative instrument normalizations. We left 
the relative XIS-BI/XIS-FI normalization free, but best-fit values were always extremely close 
to 1.00. The PIN/XIS-FI normalization was kept fixed at 1.13, a value derived using Suzaku 
observations of the Crab (Ishida et al. 2006); the uncertainty on the PIN/XIS-FI normalization 
is also discussed in §3.1. All errors on one interesting parameter correspond to Ax^ = 2.71 
(with the XIS-BI/XIS-FI normalization left free). All fits were done using xsPEC v.11.3.2. The 
abundances of Lodders (2003) were used. A neutral Galactic column of 2.94 x 10^° cm~^ was 
included (Dickey, Lockman 1990). 

3.1. Preliminary Broadband Fits 

The X-ray continuum emission of Seyfert 1 and 1.5 AGN is usually dominated by a 
power-law component thought to originate from Comptonization of soft seed photons by a hot 
corona near the central black hole (e.g., Shapiro et al. 1976; Sunyaev, Titarchuk 1980; Haardt et 
al. 1994). A simple power-law (henceforth denoted the "primary power-law") over 0.3-76 keV 
yielded a very poor fit, with x^/dof (degrees of freedom) = 49708/1449. As shown in Figure 
2(a), residuals strongly indicated the need to include absorption to account for continuum 
curvature ^3-4 keV. The narrow 6.4 keV Fe Ka line was also obvious. 

We therefore added an absorbing column to the power-law, assuming systemic redshift 

^ The HIGHECUT model is of the form I{E) = exp{{Ec - E)/Ef) for E > Ec, and I{E) = 1 for < Ec, where 
E is the photon energy, E^ is the cutoff energy, and E^ is the e-folding energy. 
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and initially assuming a covering fraction of unity (this component is henceforth denoted the 
"primary absorber"). We used xstar tables that assumed an underlying continuum with a 
photon index of F = 2.0 to model absorption in this paper. We also added Gaussians to model 
Fe Ka and Kf3 emission lines. The Fe Kf3 energy centroid was kept fixed at 7.06 keV; the 
normalization was kept fixed at 0.13 times that of the Ka line. With these changes to the 
model, x^/dof fell to 8312/1444. The best-fit ionization parameter and column density were 
log(0 = 2.0±0.1 erg cm s~^ and 7.1±0.2 x 10^^ cm"^, respectively = 4nFion/n; F^on is the 
1-1000 Ryd ionizing continuum flux; n is the density of the reflecting material). However, as 
shown in Figure 2(b), this model did not accurately describe the broadband emission. 

To attempt to model the soft emission, we added a second power-law (the "soft power- 
law"), with photon index F tied to the primary power-law; /dof fell to 3039/1443. As shown 
in Figure 2(c), residuals suggested the presence of soft X-ray emission lines, e.g., near 0.56 keV, 
which is likely due to O VII. 

We added 13 Gaussians to our model, widths a were fixed at 0.5 eV. Energy centroids 
for the lower signal-to-noise ratio lines were kept fixed at lab-frame energies. Table 1 shows the 
results for the lines in our best-fit baseline model (see below). Data/model residuals for these 
emission lines are shown in Figure 3. Removing lines one at a time from the final fit suggested 
that it was significant at > 99.0% confidence in an F-test to include each line in the fit. However 
we caution that the resulting F-test significance levels were likely upper limits in a few cases, 
since there could be some blending of Gaussian profiles when two lines' energy centroids are 
relatively close together. In addition, changing the order in which lines were added might also 
affect the derived significance of any one line.^ We identify these lines as originating in H-like 
C, N, O, Ne and Mg, and He-like N, O, Ne and Mg. We also report radiative recombination 
continua (RRC) of O VII, O VIII and Ne IX (and possibly C VI, blended with the 0.500 
keV N VII emission line). The lines are likely due to photo-ionization. Most of these lines 
have been reported previously; we refer the reader to Turner et al. (2003) for results using 
the XMM-Newton-KQ'&. The strongest line detected in both the RGS spectrum and the Suzaku 
spectrum is due to O VII; line intensities of N VI, N VII/C VI RRC blend, O VII, and O VIII 
as measured by Suzaku were roughly consistent with the RGS measurements. We also included 
a line near 0.83 keV for Fe L XVII, to model any contribution from collisionally-ionized gas. 
Turner et al. (2003) included a Mg XI recombination edge component. However, we eschewed 
it in favor of a Mg XII line at 1.47 keV; the higher effective area of the XIS compared to 
that of the RGS at this energy allowed us to determine that a line gave a better fit than a 
recombination edge. Including these emission lines in the model resulted in x^/dof falling to 
2347/1424. Figure 2(d) shows the data/model residuals after the soft X-ray lines are modeled. 

Residuals in the PIN band, peaking near 20-30 keV, signaled the need to include a 

^ For further discussion regarding justification of including multiple lines in complex spectral fits, we refer the 
reader to Pounds & Vaughan (2006), who demonstrate an application of a Bayesian analysis technique. 
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Compton reflection component. We modeled this by adding a pexrav component (Magdziarz, 
Zdziarski 1995), assuming solar abundances and an input photon index tied to that of the 
primary power-law. We fixed the cut-off energy at 400 keV. This choice was somewhat arbitrary, 
but the cutoff energy is not well constrained; in the baseline model below, we found a lower 
limit of 120 keV. We initially fixed the inclination at 30°, as per NGC 3516's classification as 
a Seyfert 1.5, and we initially assumed that the reflector is subject to the same absorption as 
the primary power-law. In the best flt-model, x^/dof fell to 1895/1423, and the value of the 
reflection fraction R (deflned as f2/27r, where Vt is the solid angle subtended by the reflector) 
was 2.8lo;2. The photon index F was 2.142;'^q;oo7- As shown in Figure 2(e), the residuals in the 
PIN band were thus corrected. 

The uncertainty on the relative PIN/XIS cross- normalization is about 3% (Ishida et 
al. 2006); the subsequent effect on R is smaller than that associated with the uncertainty of 
the PIN NXB. We modifled the intensity of the PIN NXB by ±2%, which is the la level of 
the current reproducibility of the PIN NXB. R changed by ±0.2, which is smaller than the 
statistical error on R. Errors quoted on R for the remainder of this paper are the statistical 
errors only; readers should bear in mind the additional systematic uncertainty associated with 
the NXB. 

There remained residuals in the 5-6 keV band. We flrst discuss modeling these residuals 
using relativistic diskline components, and later (§3.2) we will discuss if a partial covering 
component could provide as good a flt. We added two relativistic diskline components for Fe 
Ka and Fe K/3 emission, using a Kerr (maximally rotating) black hole line proflle (Laor 1991). 
Formally, the K/3 diskline is not required (omission of this component does not change /dof 
signiflcantly) , but we include it for completeness. The normalization of the K/3 diskline was 
flxed at 0.13 times that of the Ka diskline. The Fe Ka line energy was constrained to lie within 
rest-frame energies of 6.40 and 6.96 keV; the K/3 line energy was flxed at 7.06 keV. All other 
parameters were kept tied between the Ka and K/3 components. The emissivity index (3 (used 
when quantifying the radial emissivity per unit area as a power-law, r~^) was flxed at 3.0. 
Initially, the disk inclination i was flxed at 30°. The outer radius -Rout was kept flxed at 400 
i?g (1 i?g = GMbh/c^)- The inner radius i?in was left free. With the best-flt model, /dof 
fell to 1643/1420 (Ax^ = -252), and the data/model residuals near 5-6 keV were improved 
considerably; see Figure 2(f). 

Refltting with the diskline inclination i as a free parameter yielded a signiflcant improve- 
ment in the flt: x^ldof fell to 1521/1419 (Ax^ = -122) for i < 23°, signiflcant at >99.99% 
confldence in an F-test. In this and all subsequent best-flt models, we flxed the inclination of 
the Compton reflection component to match that of the diskline as opposed to leaving it flxed 
at 30°, though in practice this usually had negligible impact on the flt. We refer to this model 
as the "1-absorber + Compton reflection + diskline" model. Additional model parameters are 
listed in Table 2. 
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Visual inspection of the data/model residuals in the Fe K bandpass suggested an addi- 
tional dip near 6.9 keV, at the rest-frame energy for Fe XXVI; see Figure 4(a). H- and He-like 
Fe K absorption features might be expected, given the detection of a high-ionization absorber 
by Turner et al. (2005). Adding a Gaussian at 6.96 keV, fell by 7.9 (for one less dof), 
significant at 99.4% in an F-test. The absolute values of the intensity and equivalent width 
EW relative to the primary power-law were 3±2 x 10"^ ph cm~^ s~^ and 9±5 eV. Adding a 
narrow (cr = 0.5 eV) Gaussian with energy centroid fixed at 6.70 keV yielded an improvement 
in fit of Ax^ = -6.3, significant at 98.4% in an F-test. The absolute values of the intensity and 
equivalent width EW were 4±2 x 10~^ ph cm~^ s~^ and 9±4 eV, respectively. 

The Chandra-HETGS spectrum yielded narrow absorption features in the 1-3 keV band 
due to highly-ionized Mg, Si, and S, but we do not significantly detect any narrow absorption 
features at those energies (and we would likely not expect to, given the XIS resolution). A small 
dip in the spectrum near 2.3 keV is close to the expected energy for S XII. However, xstar 
models demonstrate that S XV absorption, though not significantly detected here, is always 
stronger than S XII. This feature is more likely due to calibration uncertainty associated with 
an instrumental Au M edge. 

To model Fe K absorption features, we added a second absorbing component (henceforth 
denoted the "high-ionization absorber"), again using an xstar table. Based on the results of 
Turner et al. (2005), we assumed an outflow velocity of 1100 km s~^. /dof fell to 1485/1417 
(Ax^ = -36). In this model, the high-ionization absorber had a column 4.0^*^3;! x 10^^ cm~^, 
similar to the value used by Turner et al. (2005), and log((^) = 3.7^^7 erg cm s~^. We henceforth 
refer to this model as our "2-absorber" or "baseline" model. Data/model residuals are shown in 
Figures 2(g) and 4(b); the 6.9 keV residuals are reduced and there is some slight improvement to 
the >7 keV continuum. The primary power- law, with F = 1.9ll[j:[]^, was absorbed by a column 
5.5±0.2 X 10^^ cm~^ and log((^) = 0.3±0.1 erg cm s~^. Forcing the ionization parameter to a 
much lower value (e.g., log(.^) ~ -0.5) resulted in a significantly worse fit, with large residuals 
at 1.5-3.0 keV; this is likely because the lower ionization does not account for absorption edges 
due to higher ionization species of Si and S. We caution, however, that this result could be 
influenced by residual calibration effects in the XIS near 1.8 keV and 2.3 keV (instrumental Si 
K and Au M edges, respectively). Refitting the model with the diskline emissivity as a free 
parameter yielded no additional improvement to the fit. We found (3 = 3.2±0.3; we leave f3 
frozen at 3.0 in subsequent fits. The values of EW for the broad and narrow lines were 287^24 
and 123±7 eV, respectively. Figure 4(c) shows the data/model residuals when the broad lines 
were removed from the model. When the model was then re-fit, there were correlated residuals 
in the Fe K bandpass, as illustrated in Figure 4(d). In addition, the resulting high value of 
X^/dof (1651/1421; Ax^ = +166) compared to the baseline model indicated that removing the 
disklines yielded a significantly worse fit (>99.999% in an F-test).^ 

We note that when removing the PIN data, and fitting using only the XIS, the uneertainty on the broad 
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In the best-fit baseline model, the strength of the Compton reflection hump was R = 
I.TIqIs; ^ contour plot of i? as a function of T is shown in Figure 5. Also shown in Figure 5 is a 
contour plot of i? as a function of the PIN/XIS-FI normalization, which had been kept frozen 
at 1.13; in this plot, F was a free parameter. Figure 6 shows an unfolded model spectrum. 
Other model parameters are listed in Table 2 (see Table 1 for the soft X-ray emission lines). 

Finally, we re-fit the model, assuming that the Compton reflection component was not 
affected by the primary absorber. This yielded a goodness of fit nearly identical to the previous 
fit, with x^/dof = 1479/1417. All fit parameters were virtually identical to the previous fit; for 
instance, R was 1.9±0.4; we will continue to assume that the Compton reflection component is 
absorbed. 

We returned to focus on the soft X-ray emission lines, and attempted to model all emis- 
sion lines using an xstar grid for photo-ionized line emission; the grid assumed an underlying 
photon index of 2.0. We also included an extra narrow Gaussian at 0.84±0.02 keV, again 
likely associated with Fe L XVII emission from coUisionally-ionized gas. Assuming a single 
photo-ionized emission zone with abundances for N, O, Ne, Mg and Fe left as free parameters 
yielded x^/dof = 1559/1424 for values of the ionization parameter log((^) = 1.3lo;2 erg cm s~^ 
and column density 1.4l'^o;5 x 10^^ cm~^. Best-fit abundance values relative to solar were = 
1.2tJ:^, Zo = 0.5toi, Znc = l-Otoi, Zug = 2.8±0.6, and Zpe = 0.3±0.3, i.e., there is a similar 
indication of a high value of N/0 as inferred by Turner et al. (2003) from the XMM-Newton 
RGS spectrum. Adding a second zone of emission from ionized gas with a significantly different 
ionization parameter did not improve the fit. While the fit using the one-zone model is not 
poor, for the remainder of this paper, we will continue to use the multiple Gaussians to model 
the photo-ionized emission lines, as that model yielded a lower value of x^/dof. 

3.2. Additional Partial Covering Components 

We next explored the possibility of an additional, partial covering component whose 
presence could potentially be manifested as curvature in the 1-5 keV continuum. We differen- 
tiate such components from the soft power-law, which itself could indicate "leaked" emission 
as part of a partial covering scenario (see §7.1). Starting with the baseline model, we added 
a partial covering component consisting of a power-law (with photon index tied to that of the 
primary power- law) absorbed by an xstar component. We first kept the column density A^h,pc 
and ionization parameter log(^pc) tied to those of the primary absorber. In the best-fit model 
(henceforth denoted Model PCI), the added power-law had a normalization 0.32±0.03 times 
that of the primary power-law. As shown in Table 3, A^h and log(^) of the primary absorber and 
the diskline parameters were consistent with values of the baseline model. However, y^/dof 
was 1486/1416, virtually identical to the baseline model, and so there is no formal requirement 
to include the new partial covering component when the column density and ionization states 

line flux increases by 20%, and the uncertainty on F increases by 80%. 
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are tied to those for the primary absorber. 

Next, we untied A^h,pc and log((^pc) and refit (Model PC2). The best-fit parameter 
values are listed in Table 3. Compared to Model PCI, dropped by only 5.2 for 2 additional 
degrees of freedom, significant at only ~90% in an F-test. It is thus not highly significant to 
include the new partial covering component in this case either. 

Next, we addressed whether it was possible for a partial covering component to mimic 
the curvature in the ~4-6 keV continuum modeled above as a relativistically-broadened Fe 
line. Starting with the baseline model, we removed the disklines, and added a partial covering 
component consisting of a power-law (with photon index tied to that of the primary power- 
law) absorbed by gas with a relatively low value of the ionization parameter (Model PCS). The 
column density needed to be ^10^^'^ cm~^ in order to have a rollover near 2-3 keV, causing 
an apparent "peak" near ~4-6 keV (we forced it to be greater than 4 x 10^^ cm~^). In the 
best-fit model, the partial coverer had a column density of 7.8^*^3^8 x 10^^ cm^^ (error pegged at 
lower limit), and log(^pc) = -0.4l2!6 (error pegged at lower limit). The new partially-covered 
power-law had a normalization <0.22 times that of the primary power-law. However, x^/dof 
was 1677/1418 (Ax^ = +192), a much worse fit compared to the baseline model, and there 
were still large data/model residuals in the Fe K bandpass; see Figure 2(h). We conclude that 
a partial coverer cannot mimic the observed curvature of the diskline. 

3.3. Relativistic Reflection Fits 

The broad Fe K diskline component is a signature of reflection off a possibly-ionized 
disk. Suzaku's broad bandpass makes it an ideal instrument for attempting to model the 
entire reflection spectrum (broad Fe line plus hard X-ray reflection continuum plus soft X-ray 
emission) of Seyferts in a self-consistent manner. Speciflcally, we used the ionized reflection 
models of Ross, Fabian (2005) modifled by relativistic smearing. 

We flrst removed the disklines and Compton reflection components from the baseline 
model, and modeled a blurred reflector by convolving a Ross, Fabian (2005) reflection model 
with the kernel associated with a maximally rotating Kerr black hole line proflle model (Laor 
1991) using KDBLUR. The photon index of the illuminating continuum was tied to that of the 
primary power-law. Free parameters here were the emissivity index inner disk radius Ri^ 
and disk inclination i. We also included the narrow Fe K line at 6.4 keV; this model assumes 
that all of the observed Compton reflection is associated with the broadened Fe K line and 
none with the narrow Fe K line. The best-flt model had x^/dof = 1486/1416, similar to the 
best overall flt; the results suggested a fairly neutral or lowly-ionized reflector, with log((^) 
< 1.8 erg cm s~^. The best-flt value of R was 3.0±0.8; the photon index P was 2.06±0.08. 
Parameters associated with blurring were similar to what was obtained in the baseline model: 
an inclination of 26±7° and an inner disk radius < 5.5 Rg (the emissivity index was flxed at 
3.0). The absorber parameters were consistent to those obtained for the baseline model. 
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Next, we tested a model with both a blurred reflector and an unblurred reflector, to 
test the notion that both the broad and narrow lines each be associated with contributions 
to the observed Compton reflection continuum. The blurred reflector was again modeled by 
convolving a Ross, Fabian (2005) reflection model with with a relativistic diskline proflle. The 
unblurred reflector consisted of a narrow 6.4 keV line and a Compton reflection hump modeled 
using PEXRAV. The flt was very similar to the previous one and to the baseline model: "X^ /dof 
= 1485/1415. The strengths of the blurred and unblurred Compton reflection components 
were 1.8±0.7 and 1.6±0.5, respectively. All other parameters were consistent with those of 
the previous flt (e.g., the EW of the narrow Fe K line was identical to that measured in the 
baseline model). 

3.4- Narrow Fe K Line Properties 

In our baseline model, the best-flt energy centroid for the narrow Fe Ka line was 
6.398±0.004 keV, consistent with neutral Fe. The observed line width o"obs was 26^13 ^V. 
The intrinsic line width cTintr was found by subtracting in quadrature the ^^Fe calibration line 
width (J of <4 eV from the measured line width. We inferred a 99% confidence limit for two 
interesting parameters of aintr < 45 eV, which corresponds to a FWHM velocity width <4900 
km s~^. Figure 7 shows a contour plot of line intensity versus FWHM velocity width. This 
width is consistent with the results obtained by C/iandra-HETGS for the two observations in 
2001, 1290lii[J and 3630l?i^[| km s'^ (Yaqoob, Padmanabhan 2004). 

Figure 8 shows a contour plot of the broad line intensity versus the narrow line intensity, 
illustrating that the two lines are detected independently at >4cr confidence. Such a result is 
a product of the combination of the narrow response of the XIS (yielding extremely high sig- 
nal/noise in the narrow line) and Suzaku's broad bandpass. Similar results have been reported 
e.g., for the Suzaku observation of NGC 2992 (Yaqoob et al. 2007). 

Finally, we discuss limits to a Fe Ka Compton shoulder. We added a Gaussian emission 
line at 6.24 keV (rest-frame), with width tied to that of the Ka core. We found an upper 
limit to the intensity of 7 x 10~^ ph cm~^ s~^, or 13% of the Ka core intensity. This limit 
corresponds to an EW of 21 eV. 

4. Timing Analysis 

To compare the variability properties of the primary and soft power-laws, we extracted 
light curves, summed over all four XlSes, orbitally-binned, and background-subtracted, for the 
0.3-1.0 and 2-10 keV bands. The 2-10 keV light curve is plotted in Figure 9, along with the 
12 — 76 keV background-subtracted PIN light curve, binned every three orbits. The 0.3-1.0 keV 
light curve was consistent with being constant within the errors and is not plotted. We note that 
the PIN errors shown in Figure 9 are statistical only, and do not take into account systematic 
uncertainty associated with subtraction of the NXB component. For instance, Mizuno et al. 
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(2006) noted that the systematic error is roughly 6% for the 15-40 keV band over a 5760 s 
bin. We calculated the fractional variability amplitude Fvar (which quantifies the variability 
in excess of measurement noise) and its uncertainty following Vaughan et al. (2003). For the 
2-10 keV band, Fvar was 9.2 ±0.3%, with a maximum/minimum flux ratio of roughly 1.4. For 
the 12-76 keV band, no significant variability in excess of that due to measurement errors was 
detected, with an upper limit of 4.4%. The upper limit on Fvar for the 0.3-1.0 keV band was 
2.5%. 

5. Flux-resolved Spectral Fits 

We performed flux-resolved spectral fits to search for any physical connection between 
the soft and primary power-laws. Despite the limited flux range exhibited during this obser- 
vation, we attempted to determine, e.g., if the observed X-ray flux variability could be due to 
rapid variations in column density of the primary absorber, if both power-laws vary together, or 
if one power-law is constant. We split the time-averaged spectrum into periods when the 2-10 
keV flux was higher and lower than the average 2-10 keV flux of 2.31 x 10~^^ erg cm~^ s~^, as 
illustrated in Figure 9. Net exposure times for the high and low-fiux spectra for each XIS (and 
the PIN) were approximately 62.1 (46.3) and 71.3 (59.5) ks, respectively. The average 2-10 
keV fluxes were 2.47 and 2.16 x 10^^^ erg cm^^ s^^, respectively. We applied the best-fitting 
baseline model from the time-average spectrum to both spectra. All narrow Gaussian energy 
centroids and widths, jS, i and Rin for the disklines, and log(^) for the high-ionization absorber 
were kept frozen at their time-averaged values. However, all freed parameters (including F and 
R) were consistent at the 90% confidence level. For instance, we find no strong evidence that 
the column density of the high-ionization absorber varies on short timescales. A broadband 
observation of NGC 3516 spanning a larger flux range is thus needed to potentially distinguish 
determine if the two power-laws vary in concert. 

6. The Search for Red and Blue-shifted Narrow Lines 

We searched for additional narrow absorption or emission features in the Fe K bandpass, 
as seen so far in several Seyferts, including NGC 3516. For instance. Turner et al. (2002) found 
emission lines near 5.57, 6.22, 6.41, 6.53 and 6.9 keV in the November 2001 Chandra-HETGS 
and XMM-Newton EPIC spectra of NGC 3516. In the April 2001 observation, Iwasawa et al. 
(2004) claimed the presence of a transient feature with varying energy and with flux varying 
on timescales of 25 ks. One interpretation was that such features were red- and blue-shifted Fe 
K lines associated with transient "hot spot" emission on the inner accretion disk. 

We searched for such features by adding a Gaussian component to the time-averaged 
spectrum, with width frozen at 0.5 eV, sliding it over 4-9 keV in energy. In addition to the 
aforementioned absorption line at 6.96 keV, there were only two "candidate" feature with Ax^ 
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< -5.0, absorption lines near 6.0 and 6.7 keV. However, we performed Monte Carlo simulations 
to assess the statistical significance of these features (see §3.3 of Porquet et al. 2004 and §4.3.3 
of Markowitz et al. 2006 for a description of these simulations, also Gallo et al. 2005), and we 
found that the lines were consistent with photon noise. Furthermore, each candidate features 
was evident only in one XIS camera, and thus was likely not real. 

As an aside, we note that a 7.47 keV Ni Ka emission line was not detected in the 
time-average spectrum; adding a Gaussian at this energy, we found an upper limit of 3 eV. 

We also extracted time-resolved spectral slices by dividing the time-averaged spectrum 
into five slices 48 ks in duration, with each slice having an exposure time near 27 ks per XIS and 
21 ks for the PIN. We did not investigate longer time slices since they might miss short-lived hot- 
spot emission lines; shorter time slices would have yielded poorer photon statistics. Applying a 
sliding narrow Gaussian over 4-9 keV in each spectral slice revealed only 3 "candidate" emission 
or absorption lines (with Ax^ < -6). Again, however, the features were seen in only one or 
two XIS cameras, and Monte Carlo simulations showed that not a single candidate feature 
was inconsistent with photon noise at greater than 80% confidence. Analysis of the high-and 
low-fiux spectra similarly yielded no significant narrow emission or absorption lines (even at 
6.70 and 6.96 keV). Typical upper limits to the equivalent width of an emission line at 5.57 
keV (one of the energies of the transient lines in Turner et al. 2002), were ^10-15 eV in the 
time-averaged spectrum or any of the sub-spectra. 

7. Discussion 

During the late 1990's, NGC 3516 typically displayed a 2-10 keV flux of ~4-6 x 10-^^ 
erg cm~^ s~^ (e.g., Markowitz, Edelson 2004). During the 2001 XMM-Newton/ Chandra ohser- 
vations, however, the observed 2-10 keV flux was much lower: 1.6-2.3 x 10~^^ erg cm~^ s~^ 
(Turner et al. 2005). Table 4 lists the inferred absorption-corrected 2-10 keV nuclear fluxes from 
the XMM-Newton observations, as well as during the 2005 Suzaku observation. In addition. 
Figure 10 shows the unfolded observed spectra for the Suzaku XIS and the 2001 XMM-Newton 
EPIC-pn data (see Turner et al. 2005 for details regarding the XMM-Newton data). The Suzaku 
observation apparently caught the source in a similar low level of nuclear flux as the 2001 ob- 
servations. The observed 0.5-2.0 keV flux during the Suzaku observation, however, was ~2-3 
times lower than during the XMM-Newton observations, indicating that the source was still 
heavily obscured, and confirming that the complex absorption in this source cannot be ignored 
when fitting the broadband spectrum and modeling diskline components. 

7.1. Power-law Components 

The primary power-law observed in the hard X-rays is likely emission from a hot corona 
very close to the supermassive black hole, as seen in all Seyferts. The nature of the soft 
power-law component, however, is not as clear. It could represent nuclear emission scattered 



14 



off optically-thin material, e.g., in the optical/UV Narrow Line Region (NLR). In the baseline 
model, the normalization of the soft power-law relative to that of the primary power-law was 
4.2±0.4%. Assuming a covering fraction of unity, this ratio is equal to the optical depth of 
the scattering material, indicating a column density of roughly 5 x 10^^ cm~^, consistent with 
this notion, though the column density is somewhat too high to likely be associated with the 
NLR. It is interesting to note that this column density is similar to that obtained for the 
high-ionization absorber, suggesting the possibility that this absorbing component could be 
associated with a zone of scattering. Using Chandra- ACIS, George et al. (2002) found the 
extended circumnuclear gas to have a 0.4-2.0 keV flux of roughly 10~^^ erg cm~^ s~^. However, 
that emission was studied over an annular extraction region 3" to 10" (0.6-1.8 kpc), and so 
that flux value is likely a lower limit to the 0.4-2.0 keV flux that Suzaku would observe (given 
the XRTs' 2' HPD). In our baseline model, we found an unabsorbed 0.4-2.0 keV flux of 1 x 
10~^^ erg cm~^ s^^, consistent with the notion that the soft power-law is scattered emission. 

In this case, the decrease in observed 0.5-2.0 keV flux from 2001 to 2005 could potentially 
be explained by the scattered emission responding to a recent decrease in nuclear continuum 
flux. However, this scenario would require the bulk of the scattered emission to lie within 
at most a few light years of the black hole, and the nuclear flux would have had to decrease 
between 2001 and 2005 (when the source was not observed by any major X-ray mission in the 
2-10 keV band) then return to 2001 levels by the Suzaku observation. 

Alternatively, the soft power-law could be unobscured, "leaked" nuclear emission as part 
of a partial covering scenario. In this case, the primary absorber would obscure 96% of the sky 
as seen from the nuclear continuum source. The lack of significant variability in the 0.3-1.0 
keV band could argue for the soft power-law to originate in scattered emission, since we might 
expect to observed variability of the same amplitude as the 2-10 keV band only if the soft 
power-law were leaked nuclear emission. However, this is far from certain, as the 0.3-1.0 keV 
band had a low count rate and the presence of the soft emission lines in the XIS spectrum could 
contribute to dilution of intrinsic variability in the soft power-law. A broadband observation 
spanning a larger observed flux range is needed to clarify this issue. The soft power-law could 
of course represent a blend of scattered emission plus leaked nuclear emission. We therefore 
conclude that the primary absorber has a covering fraction between 96-100%. 

7.2. Complex Absorption 

We detect two zones of absorption: in addition to the primary absorber, which has 
a covering fraction of 96-100%, there is the high-ionization absorber, which is assumed here 
to have a covering fraction of unity. The high-ionization absorber is potentially the same as 
that reported by Turner et al. (2005); we find a column density A^h of 4.0^3;! x 10^^ cm~^, 
consistent with the column density of 2 x 10^^ cm~^ used by Turner et al. (2005), although we 
use a slightly higher ionization parameter (see Turner et al. 2007). Previous studies of NGC 
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3516, such as Netzer et al. (2002), have discussed in detail the UV absorber, responsible for 
H Lja, C IV and N V absorption features in Hubble Space Telescope spectra (Kraemer et 
al. 2002). In the X-ray band, discrete features associated with Mg VII-IX and Si VII-IX 
are expected from this component, but with the CCD resolution and with calibration- related 
artifacts near 1.7-1.8 keV in the XIS, such features are not detected by Suzaku. 

Suzaku has found the primary absorber of the hard X-ray continuum to be lowly-ionized 
(log(^) = 0.3±0.1 erg cm s~^), with a column density A^h of 5.5±0.2 x 10^^ cm~^. It is possible 
that it is the same absorber that Turner et al. (2005) designated as the "heavy" partial-covering 
absorber, though we use a somewhat lower ionization parameter (see Turner et al. 2007). A 
new observation of NGC 3516 with XMM-Newton in 2006 October showed that the source had 
returned to similar >6 keV brightness and similar obscuration levels (A^h ~ 2 x 10^^ cm~^; 
covering fraction ~ 45%) as during the 2001 observations (Turner et al. 2007). Long-term 
changes in the covering fraction of the heavy absorber could explain the bulk of the spectral 
variability changes between the 2001, 2005 and 2006 observations. In this case, the column 
density has decreased by a factor of 4.5, while the covering fraction has increased from ~40- 
60% to 96-100%, from 2001 to 2005, and subsequently returned to approximately the same 
levels as in 2001 within the next 12 months. However, because we have not actually observed 
the entire eclipse associated with a specific, discrete blob of absorbing gas traversing the line 
of sight, it is not clear whether the covering fractions derived are associated with single, large 
blobs partially blocking the line of sight to the X-ray continuum source, or if the absorber 
consists of numerous, discrete blobs or has a filamentary or patchy structure. On the other 
hand, given the gaps between the 2001 and 2006 XMM-Newton and 2005 Suzaku observations, 
it is certainly plausible that these observations could have caught independent, discrete blobs 
or filaments with differing column densities and differing physical sizes and/or radial distances 
lying on the line of sight. 

To estimate the distance r between the central black hole and the absorbing gas, we 
can use a definition of the ionization parameter ^ = i>i_ioooRyd/(^^^); where n is the number 
density. Li_ioooRyd is the 1-1000 Ryd illuminating continuum luminosity, and the value of the 
ionization parameter is taken to be the value of 2 erg cm s~^ measured above. We estimate the 
maximum possible distance to the material by assuming that the thickness Ar must be less than 
the distance r. The column density A^h = nAr, yielding the upper limit r < Li-ioooRyd/(A^HO- 
We estimate the 1-1000 Ryd flux from the baseline model to be 9.9 x 10"^^ erg cm~^ s~^, 
which corresponds to Li_ioooRyd = 1-7 x 10"*^ erg s~^ (assuming Ho = 70 km s"^ Mpc"^ and 
Ao = 0.73). r is thus <2 x 10^° cm (180 light-years), a very loose upper limit encompassing 
both distances associated with the BLR (~10 light-days; Peterson et al. 2004) and a possible 
cold molecular torus at a 1 pc radius. Variability in the absorber properties between 2001-2005 
and 2005-2006 thus imply a radial distance of at most a few light years in the case of clouds 
traversing the line of sight to the nucleus. In addition, in the case of a partial covering scenario. 
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it is plausible that the absorber's size could be of the same order as that of the X-ray continuum 
source. The absorbing material could thus be e.g., associated with the base of an outflow or 
from dense clouds associated with magnetohydrodynamic disk turbulence (e.g., Emmering et 
al. 1992). 

NGC 3516's transition from an unobscured source to a moderately-obscured source in a 
4 year span presents a challenge to standard Seyfert 1/2 unification schemes. If the obscuration 
in NGC 3516 is associated with an equatorial molecular torus usually invoked in Seyfert 1/2 
unification schemes, then it is possible that during the Suzaku observation, the inner edge of 
the torus could have intersected the line of sight, but given NGC 3516's classification as a 
Seyfert 1, this could only occur if the torus opening angle were extremely small and the torus 
were not azimuthally symmetric. Alternatively, variations in column density and/or covering 
fraction could be due to fine structure in large-distance (tens of pc), non-equatorial filaments 
that traverse the line of sight (e.g., Malkan et al. 1998). 

7.3. Fe K Emission Components and Compton Reflection 

We have deconvolved the broadband emitting components, and determined that 1) the 
existence of the broad Fe line is robust in that it was required in all models for an adequate fit, 
and 2) a partial covering component could not mimic the curvature associated with a relativistic 
broad line. We note, for instance, that if we remove the diskline components from the baseline 
model and refit, not only is the fit worse (x^ increases by over 170), but the value of R becomes 
~ 3.2. This value is incompatible with the observed EW of the narrow line unless the Fe 
abundance is extremely sub-solar (^0.3; see below). The best-fit disk inclination was typically 
^25°. The inner radius was typically ^5i?g. The line energy was seen to be consistent with 
neutral to mildly- ionized Fe (up to Fe ~ XX; Kallman et al. 2004). The equivalent width with 
respect the primary continuum was 287^24 eV, consistent with the value of 43l!t}72 obtained 
by Turner et al. (2005) for the April 2001 XMM-Newton observation, where the spectrum could 
be fit with a diskline component in addition to the complex absorbing components. 

The line energy of the narrow Fe Ka line was also consistent with emission from neutral 
Fe. The intensity of the narrow line during the Suzaku observation is roughly 40% higher than 
that measured during the 2001 XMM-Newton and Chandra observations (Turner et al. 2002), 
possibly indicating that a substantial fraction of the Fe K line photons originate in a region 
^5 It.-yr. in size. We measured a FWHM velocity line width for the narrow Fe Ka line of < 
4900 km s~^ (99% confidence level for two interesting parameters). This velocity does not rule 
out an origin in the BLR; Peterson et al. (2004) reported FWHM velocities for the Ha and 
H/3 lines of 4770±893 and 3353±310 km s~^, respectively. However, we also cannot exclude a 
contribution from an origin in the putative molecular torus; there could potentially be a very 
narrow line component with FWHM velocity ~ a few hundred km s~^, but the XIS would be 
unable to separate it from the relatively broader line component. 
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It is possible that the same material that absorbs the hard X-rays along the line of sight 
is responsible for producing the narrow Fe line. The material producing the Fe line cannot have 
a column substantially less than 10"^^ cm~^ or else there would be insufficient optical depth 
to produce a prominent Fe K line. The primary absorber, with its column density of 5.5 x 
10^^ cm~^ and low ionization state, is thus a plausible candidate for the narrow line origin. As 
an estimate of the Fe Ka equivalent width expected in this case, we can assume an origin in 
optically-thin gas which completely surrounds a single X-ray continuum source and is uniform 
in column density, and use the following equation: 

EWcslc = fcUjfKaA ^ — r (1) 

^ I -Celine j 

Emission is assumed to be isotropic. Here, fc is the covering fraction, initially assumed 
to be 1.0. u is the fluorescent yield, 0.34 (Kallman et al. 2004). /ko is the fraction of photons 
that go into the Ka line as opposed to the K/3 line; this is 0.89 for Fe I. A is the number 
abundance relative to hydrogen. We assumed solar abundances, using Lodders (2003). P{E) is 
the spectrum of the illuminating continuum at energy E; E^j^c is the Ka emission line energy. 
(Tph(-E) is the photo-ionization cross section assuming absorption by K-shell electrons only 
(Veigele 1973^). 

For A'h = 5.5 x 10^^ cm~^, EWcaic = 29 eV, substantially lower than the observed 
EW of 123±7 eV. We conclude that it is possible for the primary absorber to contribute to 
the observed line EW, but there is also likely a contribution from some other (non-continuum 
absorbing) material lying out of the line of sight, likely with column densities 10~^'^ cm^^ (e.g., 
Matt et al. 2002). For instance, if the putative cold molecular torus does not intersect the line 
of sight, it could contribute to the observed EW. The 13% upper limit to ratio of the Compton 
shoulder/ narrow Fe Ka core intensity was a 90% confidence limit only, and does not exclude 
at high confidence the possibility of Compton-thick material out of the line of sight. 

An additional possibility is that the material emitting the bulk of the line photons 
could be responding to a continuum fiux that was higher in the past. For instance, if the 
putative molecular torus is located ~ a pc or so from the black hole, the torus will yield a line 
EW corresponding to the continuum fiux averaged over the past few years. This situation is 
plausible for NGC 3516, as the 2-10 keV fiux of NGC 3516 during ^1998-2001 (Markowitz, 
Edelson 2004) was a factor of ~1.5-2 times brighter than during 2005. 

We now turn our attention to properties of the Compton refiection continuum. Suzaku 
has observed other Seyferts to display reduced levels of variability in the PIN band compared 
to the 2-10 keV band, e.g., in MCG-6-30-15 (Miniutti et al. 2007). This behavior is thought to 
be caused by the presence of the relatively non- varying Compton refiection hump, which dilutes 
the observed >10 keV variability of the power-law component. Gravitational light-bending in 
the region of strong gravity has been invoked to explain the relative constancy of the refiection 

^ http://www.pa.uky.cdu/^vcrncr/photo.html 
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spectrum (Compton hump and Fe K diskline component) despite large variations in the coronal 
power-law flux in MCG-6-30-15, for instance (Miniutti et al. 2007). In the case of NGC 3516, 
the observed fractional variability amplitudes for the 2-10 and 12-76 keV bands were Fvar,2-io 
= 9.2±0.3% and Fvar, 12-76 < 4.4%, respectively. These measurements allow us to rule out the 
possibility that the Compton hump varies in concert with the power-law, since the variability 
amplitudes would be consistent in that case. The primary power-law and Compton hump 
contribute 44% and 56%, respectively, of the total 12-76 keV flux. In the case of a constant 
Compton hump and variable power-law, Fvar, 12-76 would then be equal to Fvar,2-io / 2.25, or 
roughly 4.1%. This is roughly consistent with the observed upper limit on F^ar, 12-76) suggesting 
that the reflection component varies less strongly than the primary power-law over the course 
of the observation. To verify this, however, we would need to observe NGC 3516 over a larger 
X-ray flux range than in the current Suzaku observation to potentially observe any significant 
variability in the PIN band. 

Finally, we discuss the origin of the material that gives rise to the observed Compton 
reflection hump. The primary and high-ionization absorbers lack the necessary column density, 
and are excluded. We next consider an origin in the same material that yields either the 
broad or narrow Fe lines. George, Fabian (1991) calculated that R = 1 corresponds to an 
observed Fe Ka line EW (relative to a primary continuum with a photon index of 1.9) of 140 
eV for neutral Fe, assuming an inclination angle of 25°. However, George, Fabian (1991) used 
the elemental abundances of Morrison, McCammon (1983), where the Fe number abundance 
relative to hydrogen was A-pc = 3.3 x 10~^. More recent papers have slightly lower values of A-pe, 
2.7-3.0 X 10^^ (Lodders 2003; Wilms et al. 2000). The expected equivalent width corresponding 
to i? = 1 is thus 115-125 eV. In our baseline model, we found a best-fit value oi R = 1.7, which 
corresponds to an expected line EW (relative to the primary continuum) of 200-215 eV, a value 
in between the observed EWs of the broad line (287 eV in the baseline model) and the narrow 
line (123 eV). It is thus not clear from this measurement alone whether the total Compton 
reflection continuum is associated with the broad line (disk), narrow line (a distant origin), or 
both. That is, while is it a possibility that at least some portion of the Compton reflection 
component is associated with the broad Fe K component, we cannot exclude the possibility that 
the narrow line contributes as well and that there is reflection off cold, distant material. For 
example, in §3.3, we demonstrated that a model wherein there existed both blurred reflection 
from an ionized disk plus reflection from cold, distant material, such as the molecular torus, 
provided a good fit to the data. In addition, we demonstrated in this section that the observed 
EW of the narrow Fe K line means we cannot rule out a contribution to the narrow Fe line, 
and to the reflection continuum as well, from Compton-thick material out of the line of sight. 
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8. Summary of Main Results 

We have reported on a 150 ksec observation of NGC 3516 obtained with the Suzaku 
observatory in October 2005. The good exposure times after screening were 135 ks for each of 
the XIS cameras and 106 ks for the HXD-PIN. 

Our best-fit broadband model included a primary power-law with photon index F = 
1.904±0.025 in our baseline model, absorbed by a column of material with A^h = 5.5±0.2 x 
10^^ cm~^ and with log(.^) = 0.3 erg cm s~^. We modeled the soft band continuum emission 
using a power-law component which could represent nuclear emission off optically-thin material, 
unobscured "leaked" nuclear emission, or a blend of both. The hard X-ray absorber could thus 
be a partial coverer, with a covering fraction >96%, or it could obscure the X-ray continuum 
source completely. If this absorber is the same "heavy" absorption component reported by 
Turner et al. (2005) in the 2001 XMM-Newton observations, then between 2001 and 2005 the 
column density of this absorber decreased by a factor of 4.5, while the covering fraction increased 
substantially, leading to an observed 0.5-2.0 keV flux a factor of 2-3 lower in 2005 than in 2001. 
Subsequently, by the 2006 October XMM-Newton observation, the covering fraction returned 
to approximately the same level observed in 2001. One possibility for the variations in the 
properties of the obscuring material between the 2001, 2005, and 2006 observations could be 
the presence of discrete clouds or fllaments within a few light years of the black hole traversing 
the line of sight; the equatorial molecular torus invoked in Seyfert uniflcation schemes is likely 
not directly responsible. We also modeled a highly-ionized absorber with a column density A^h 
of 4.0l3;i X 10^^ cm~^, ionization parameter log((^) = S.T^q'^ erg cm s~^, assumed to have a 
covering fraction of unity. 

Our baseline model also included a dozen narrow K-shell emission lines originating in 
He-like N, O, Ne and Mg, H-like C, N, O, Ne and Mg and three RRC features, consistent 
with an origin in photo-ionized material. However, we cannot exclude a contribution from 
collisionally-ionized material, as suggested by the presence of an Fe L-shell XVII line near 0.83 
keV. 

The broad Fe Ka line has been robustly detected: we can distinguish between the 
curvature in the observed continuum due to a partial coverer and that due to a broad diskline; 
we conclude that for this observation of NGC 3516, a diskline component is required and that 
neither a cold nor ionized partial coverer can mimic the continuum curvature associated with 
the diskline component. The broad and narrow lines are decoupled (detected independently) at 
high signiflcance, thanks to the narrow response of the XIS and the subsequent high signal/noise 
ratio in the narrow line. In our best-flt model, we flnd the Compton reflection strength to be 
the value of i? = l.lto.l The narrow Fe Ka line, meanwhile, has a FWHM velocity width of 
<4900 km s~^ (99% confldence level for two interesting parameters), consistent with an origin 
in material with the same velocity as NGC 3516's BLR, though a contribution from material 
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with lower velocity widths cannot be excluded. It is possible that the primary hard X-ray 
absorber may also be responsible for emitting the narrow Fe Ka line, though there may be a 
contribution from material lying out of the line of sight, such the putative molecular torus. 

The authors gratefully acknowledge the dedication and hard work of the Suzaku hardware 
teams and operations staff for making this observation possible and for assistance with data 
calibration and analysis. This research has made use of HEASARC online services, supported 
by NASA/GSFC. This research has also made use of the NASA/IPAC Extragalactic Database, 
operated by JPL/California Institute of Technology, under contract with NASA. 
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Table 1. Soft X-ray Emission Lines 
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Results are for our best-fit baseline model. 
^ denotes a fixed parameter. 
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Table 2. Best-fit parameters for the 1-absorber model and the baseline (2-absorber) model 
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component (modeled with pexrav) was set at 20°. 
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Table 3. Best-fit parameters for models with partial covering components 
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Units of power-law normalization are ph cm^^ g-i keV^^ at 1 keV. 
^ parameters tied to those of the primary absorber. 
The inclination angle of the Compton reflection component (modeled with pexrav) was set at 25°. 

Table 4. Comparison to the 2001 XMM-Newton observations 



Suzaku XMM-Newton XMM-Newton 

Oct. 2005 Apr. 2001 Nov. 2001 



Absorption-corrected 2-10 keV flux (erg cm^^ s^^) 
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Results for the XMM-Newton spectra were taken from Turner et al. (2005). 
Results for Suzaku are from the best-fit baseline model. 
^ denotes a fixed parameter. 
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Fig. 1. HXD-PIN spectrum. The upper panel shows the net source spectrum {gray points), the back- 
ground [lower black points), and the total (source + background) spectrum {upper black points). The PIN 
spectra have been binned such that the net spectrum has a minimum signal-to-noise ratio of 8cr per bin. 
The lower panel shows the ratio of the net source spectrum to the total spectrum. 
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Fig. 2. Data/model residuals for selected broadband models. Panel a) shows residuals to a simple pow- 
er-law. Panel b) shows the residuals after the primary absorber and narrow Fe lines are included in the 
model. Panel c) shows residuals after the soft power-law is included. In panel d), the soft emission lines 
have been modeled. In panel e), a Compton reflection component has been added (using PEXRAV in 
XSPEC). In panel f), diskline components have been added. In panel g), the high- ionization absorber 
has been added to yield our baseline model. Panel h) shows the results when the broad Fe lines are 
removed and substituted with an additional partial covering component. Black points <12 keV denote 
the (co-added) XIS-FI spectrum. Gray points denote the XIS-BI spectrum. Black points >12 keV denote 
the PIN data. Rest-frame energies are shown. All data have been rebinned with a binning factor of 5, 
though XIS data >7.2 keV have been rebinned by a factor of 25 for clarity. 
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Fig. 3. Ratio of the soft band data to a simple power-law, showing the prominent emission lines. For 
clarity, only the XIS BI spectrum is shown. Rest-frame energies are shown. 
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Fig. 4. Data/model residuals to various models are shown, focusing on the Fe K bandpass. Data have 
been rebinned by a factor of 3. Panel a) shows the residuals when the primary absorber is the only zone 
of absorption modeled. Panel b) shows the residuals to the baseline model, with two zones of absorption; 
note that residuals near 6.9 keV arc now improved. In panel c), the broad Fe lines have been removed 
from the baseline model. When the model is then refit, correlated data/model residuals appear in the Fe 
K bandpass, as illustrated in panel d). 



31 



1.8 1.85 1.9 1.95 1 1.1 1.2 

r PIN/XIS-FI Norm. 



Fig. 5. Left: Contour plot of the Compton reflection fraction R versus the photon index of the primary 
power-law. The black contours show results when the PIN/XIS-FI normalization was kept fixed at 1.13. 
The gray contours show results when the PIN/XIS-FI normalization was left free. Right: Contour plot of 
R versus the PIN/XIS-FI normalization; T was a free parameter. In both plots, solid, dashed, dot-dashed, 
and dotted lines denote 68%, 95.4%, 99.73%, and 99.99% confidence levels, respectively. 
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Fig. 6. Unfolded broadband spectrum for the baseline model. The thick black line is the total spectrum. 
Thin black lines denote the primary, absorbed power-law, the soft power-law, and the reflection component 
peaking at 20 keV. Gray dashed lines denote the narrow Fe lines. Gray solid lines in the Fe K bandpass 
are the diskline components. Gray solid lines in the soft band denote the soft emission lines. 
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Fig. 7. Contour plot showing the intensity of the narrow Ka line versus the FWHM velocity width. 68%, 
95.4%, 99.73% and 99.99% confidence levels for two interesting parameters are plotted (solid, dashed, 
dot-dashed, and dotted lines, respectively). 
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Fig. 8. Contour plot showing the intensity of the narrow Ka hne versus that of the broad Ka Hne. 68%, 
95.4%, 99.73% and 99.99% confidence levels for two interesting parameters are plotted (solid, dashed, 
dot-dashed, and dotted lines, respectively). Note that the 99.99% confidence level contours do not hit 0, 
demonstrating that the lines are decoupled at the 4fT level. 



35 



T ' ^ 

2-lOkeV 




J I I I I I I I I I I I L 



100 200 

Time (ksec after 2005 Oct 12, 14:33 UT) 

Fig. 9. Orbitally binned, background-subtracted 2-10 keV XIS light curve, summed over all four XIS 
cameras and background subtracted (top), and background-subtracted 12-76 keV PIN light curve, binned 
to every three orbits (bottom). The error bars show only the statistical error of the count rates; systematic 
uncertainty associated with subtraction of the non- X-ray background in the PIN is roughly 0.01 ct s~^. 
The mean count rates are 2.41 ct s~^ for 2-10 keV and 0.16 ct s~^ for 12-76 kcV. 
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Fig. 10. Unfolded spectra for the Suzaku and 2001 XMM-Newton observations. The black points denote 
the Suzaku XIS-FI and BI spectra, rebinned by a factor of 3. The light and dark gray points denote the 
XMM-Newton EPIC-pn spectra for April 2001 and November 2001, respectively; data are rebinned by a 
factor of 6. 
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